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Abstract 

Using a simple deterministic model for the Internet graph we show that the eigen- 
value power law distribution for its adjacency matrix is a direct consequence of the 
degree distribution and that the graph must contain many star subgraphs. 
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1 Introduction 

Recent research shows that most communications networks, like the world 
wide web, the Internet, telephone networks, transportation systems (including 
the power distribution network), and biological and social networks, belong to 
a class of networks known as small-world scale-free networks. These networks 
exhibit both strong local clustering (nodes have many mutual neighbors) and 
a small diameter (maximum distance between any two nodes) [13]. Another 
important common characteristic is that the number of links attached to the 
nodes usually obeys a power law distribution (is scale-free) as was observed 
first in empirical studies from the Faloutsos' brothers [10] and Barabasi and 
Albert [1]. In [10], the authors present data that show that the eigenvalue dis- 
tribution of the adjacency matrix of the Internet graph also follows a power 
law. A qualitative proof that this power-law is a consequence from the distri- 
bution of degrees of the network, which should contain many star subgraphs of 

* Corresponding author. Tel. -^34-93-413- 7072, Fax +34-93-413-7007. 
Email address: comellas@mat.upc.es (Francesc Comellas). 



Preprint submitted to Elsevier Science 



16 Deember 2004 




Fig. 1. (a) A simple model for the Internet graph: many degree one nodes are 
attached to a backbone graph B to form graph G. (b) Deleting the backbone graph 
B results in a non connected graph constituted by the union of star graphs. 

different sizes, was presented by the authors in [3]. Mihail and Papadimitriou 
in [11] provided an analytical proof by using the theory of random graphs. 
Fabrikant, Koutsoupias and Papadimitriou [9] provide an explanation of the 
Internet power laws based on complex multicriterion optimization. 

Many approaches to understanding small-world scale-free networks are based 
on stochastic models and computer simulations, however the use of determinis- 
tic models (although they do not allow to capture the full complexity of a real 
life network) do help in the understanding of their behavior and permit the di- 
rect determination of relevant parameters for the modeled networks. In partic- 
ular, it is possible to construct small-world deterministic graphs with different 
degree distributions matching the distribution of real networks, see [5,2,7,4]. 
We introduce in this paper a simple deterministic model, a toy model, to show 
analytically, that the observed eigenvalue power law of the Internet may be a 
direct consequence of the degree distribution of a star based structure (whose 
power law can be justified considering, for example, preferential attachment [1] 
or duplication [6] models). 



2 A star graphs based deterministic model 

Restricting our study to the Internet graph (at the router level), our model 
considers this graph as a scale-free network made by the union of a number of 
star graphs with different orders and connected through the root vertices by 
a relatively small backbone graph B, see Figure 1. This model is an approx- 
imation to the real Internet graph, as described in appendix A of [10], if we 
consider a degree distribution matching the observed results: 

Each Hi + 1-star graph (with a root node and Ui other nodes) will contribute 
with a vertex of degree Ui and Ui vertices of degree 1. The eigenvalues of this 
star graph are and 0. 

If we consider a set of star graphs whose root vertices have rank degrees 
distributed according to a power law P, the graph union of these graphs will 
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have eigenvalues with ranks distributed according to |. We note that, as we are 
considering a determinitic exact model, to relate the exponent of the discrete 
degree distribution to the (3 exponent of a continuous degree distribution for 
a random scale-free graph, a cumulative distribution should be considered, 
Pcum{k) = J2k'>k^{k')/N ~ A;^~^, where k and k' are values of the discrete 
degree spectrum, n{k') is the number of vertices of degree k' and is the 
order of the graph. 

It is possible now to give some insight on the eigenvalues of the global graph 
G obtained by joining the star graphs by using a backbone graph B. Since the 
number of vertices of B is small with respect to the total number of vertices of 
all the star graphs, the values of the spectrum will be very close to the original 
graph of stars and will therefore follow a power law. 

This last result is related to the interlacing theorem [8] which states that if 
r is a graph with spectrum Ai > A2 > ... > Ap, fi is a vertex of F, and 
A^i > > • • • > fJ'p-i is the spectrum of F \ Wi (the graph resulting from 
F after deleting vi and its associated edges) then the spectrum of F \ i;i is 
interlaced with the spectrum of F, i.e. Ai > /ii > A2 > /i2 > • • • > yWp-i > Ap. 

On the other hand, notice that this simple model ensures that if the power-law 
exponent for the degrees is /?, the corresponding exponent for graph G will be 
approximately (5/ 2, matching the results shown in [10]. 

In the next section we compute analitically the eigenvalues of G and verify 
that they are close to those of its star subgraphs. 



3 Analytical computation of the eigenvalues of G 



Consider s star graphs Si, S2, ■ ■ ■ , Sg with, respectively, ni + 1, ?t,2 + 1, . . . , n^ + l 
vertices. The backbone graph B joins the central vertices of each star. We call 
A the adjacency matrix of this graph B. This is equivalent to start with a 
graph B with s vertices and add to each vertex, respectively, rii, . . . ,ns new 
vertices of degree 1. 



To compute the spectrum of the global graph G {B plus the vertices of degree 
b O \y \y 



one) we consider ~ ^ which gives the system to solve 



Ax + b^y = Xx 
bx = Xy 



(1) 
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where = diag(6i, . . . , 6s) with hi = (1, . . . , 1), ones, i.e. fe-^ is a sx (J2i=i ^i) 
matrix and O is the null matrix with dimensions (J2i=ii^i) ^ (Z]i=i'^i)- 

If A 7^ 0, from (1) we obtain Ax + X~^b^bx = Xx where b^b is a diagonal ma- 
trix, whose elements are ni, . . . ,77,^. Thence Ax = (A/ — A~^diag(ni, . . . ,ns))x. 
Introducing = A/(A^ — rij), this last equation can be written 

dia.g{ai, . . . ,as)Ax = X. (2) 



If the backbone graph B is another star graph with s nodes (rooted in one of 



h 



T 



where b'^ = (1 11 



the star graphs), then the matrix A will be 

\bOs^, 

and Os-i is a null matrix of dimensions (s — 1) x (s — 1). Thus, the eigenvalues 
of G are the values of A wich verify (2) for this A, i.e. values which allow the 
matrix diag(ai, . . . , as)A to have 1 as eigenvalue. 

Solving det(diag(ai, . . . , as)A — xl) = we find 



s 



\s~2 ^2 



-x}" " I a;" — «! ^ = (3) 

i=2 J 



and the eigenvalues are with multiplicity s — 2, and i^JaiJ2i=2'^i with 
multiplicity 1. Hence, the final equation to solve is ai Yli=2 = 1- Its left part 
is a function of A 

/(A) = «1 U=2 = U=2 = E-=2 X^-rn A^-n. 

_ Y^s ni !li 

^i=2 (ni-ni)(A2-ni) („i_„,)(A2-ni) 

and it can be decomposed in a sum of simple fractions whose asymptotes are 
given by Aj^-t = ±y^. The points A'j where this rational function cuts the 
constant function 1 are the solutions of the equation, and they are therefore 
the eigenvalues of G we are looking for, see Figure 2. 

Each fraction ^ — can be decomposed as „ — ^ — „ , If we 

X-'-rii ^ 2^/n7(A-^/n7) 2^/n7(A+^/n7) 

consider the positive values to find the largest eigenvalues and assuming ni > 
. . . > Ug, for 2 < i < s approaching the value y/rii from the right, the function 
/(A) is negative, whereas approaching it from the left it is positive. Therefore it 
should cut the constant function 1 in the left side, so A^ < ^/nl, and A^ > ^rij+i 
as at ^rij+i there is another asymptote and the function is negative on the 
right side of y/rii+i. Therefore the eigenvalues of G verify: y/nl > X'- > y/rii+i 
and A'j = — A, where A < y/Wi — y/rii+i. At ^/nl the function is positive 
at the right of this point and negative at the left. Hence, A'^ > y/ni- 
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Fig. 2. Analytical determination of the eigenvalues of graph G, the graph obtained 
by joining several star graphs with a backbone graph. 



Note that if the values for rii are consecutive integers, the difference in the 
corresponding square roots will be very small for large rij and the eigenvalues 
for the new graph will be very near to those of the disconnected star graphs. 

Example 1 Consider the Oregon-Multi dataset of Internet at the AS level as 
in 2001 [12], and use a star graph as backbone graph B to interconnect the star 
graphs with Ui equal to I4OO, 1000, 900, 625, 575, 550, 515, 425, 415, 350, 
340, 320, 285, 250, 225, 215, 210, 180, 175 and 170. The degree distribution 
(by construction) fits the degree distribution of the Internet graph as shown 
in [12]. To calculate the spectra of this graph we solve the equation /(A) = 1 
and we obtain: 13.03, 13.22, 13. 4I, I4.48, I4.66, 14.99, 15.81, 16.87, 17.88, 
18. 43, 18.70, 20.36, 20.61, 22.68, 23.44, 23.96, 24.98, 29.97, 31.59, 37.80. 

In Figure 3 we plot the twenty highest eigenvalues in a log-log scale. The result 
shows a power law as in the real data from [12]. Observe that the slope of the 
regression model is very similar to the observational data. 




Fig. 3. Highest twenty eigenvalues for the Internet graph (Oregon-Multi), see [12], 
and the graph constructed in Example 1. In both cases they follow a power law. 
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4 Other models and conclusions 



The model considered here describes the degree and eigenvalue distributions 
of the Internet graph, however its clustering is zero as the resulting graph is 
a tree. We have performed the same analysis connecting the root vertices of 
the star graphs with a complete graph. The results are comparable to those 
obtained above (in that case the clustering of the global graph is different 
from zero but small). We note that if the backbone graph B joins complete 
graphs instead of star graphs this would lead to an eigenvalue power law with 
an exponent close to the exponent of the degree distribution (the eigenvalues 
of a complete graph with n vertices are n — 1 and —1). 

The star-based model considered here constitutes a convenient tool to study 
the Internet network since using deterministic techniques no simulation is 
needed and relevant network parameters can be directly calculated. Thus, the 
contrast of real with synthesized networks is straightforward. On the other 
hand, our study shows that it is a reasonable assumption to consider that the 
Internet graph contains many star subgraphs, as this explains the observed 
power-laws and relationship between the degree and eigenvalue distributions. 
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